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Abstract 



A new theoretical formalism for the optimal quantum control has been presented. 
The approach stems from the consideration of describing the time-dependent quantum 
systems in terms of the real physical observables, viz., the probability density p(x,t) 
and the quantum current j(x,t) which is well documented in the Bohm's hydrodynamical 
formulation of quantum mechanics. The approach has been applied for manipulating the 
vibrational motion of HBr in its ground electronic state under an external electric field. 
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1 Introduction 



Manipulating the outcome of a chemical dynamics by properly tailoring an external 
field has been one of the active field of research in the recent days [1-14]. Problems where 
the external field is an electromagnetic field have received the most attention [2-14], al- 
though other applications may arise as well. Theoretically, there are two basic paradigms 
for such control : a static control scheme [13,14] and a dynamic control scheme [2-12]. In 
the static scheme [13,14] one uses two or more cw light fields(optical coherence) and the 
superposition of two or more eigenstates(molecular coherence) to cause interference be- 
tween different plausible pathways to a final quantum state, and the outcome is controlled 
by tailoring different parameters of the optical and molecular coherences. Whereas the 
dynamic scheme [2-12] creates non-stationary states of one's choice, by optimally design- 
ing the electric field. This comes under the domain of the optimal control theory [15], a 
mathematical tool commonly used in the engineering fields. A basic difficulty in attaining 
the control designs is the computational effort called for in solving the time-dependent 
Schroedinger equation, often repeatedly in an iterative fashion over an extended spatial 
region. 

In this paper, we introduce a new formulation aiming at reducing the effort for the 
quantum optimal control(QOC). Our recent work [16,17] has shown that the Bohmian 
quantum hydrodynamics(BQH) is capable of being much more efficient than the conven- 
tional method(e.g., FFT propagation) and this should carry over to the optimal control 
task. This paper will show how the BQH can be utilized in the QOC. The formulation is 
based on the hydrodynamic description of the quantum mechanics emerging mainly from 
the work of David Bohm [18,19] where the dynamics is described by two equations, viz., 
the equation of motion for the probability density, p(r, t) and that for the quantum current, 
j(r,t) which are defined as p(r,t) = **(r,*)*(r,*) and j(r,t) = ±£/m[tf*Vtf - tfVtf*], 
* being the complex wave function in the time dependent Schroedinger equation(TDSE) 
and Im refers to the imaginary value. Thus one by-passes the explicit use of the time 
dependent Schroedinger equation(TDSE) and hence the typically oscillatory nature of 
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the complex wave finction. This seems benificiary at the first place because (i) one deals 
with the real quantum mechanical variables, and (ii) density and quantum current posses 
a smooth spatial variation as opposed to the wave function. Recent illustrations [16,17] 
have demonstrated the smooth spatial and temporal nature of the variables and the abil- 
ity to discretize them on a relatively small number of grid points. In pursuing the paper 
we maintain the following layout. In section 2 we give a brief account of the BQH. In 
section 3 we provide the QOC formulation based on the BQH. In section 4 we apply the 
method for manipulating the vibrational motion of HBr molecule in its ground electronic 
state. Section 5 concludes the paper. 



2 Bhomian Quantum Hydrodynamics 

Despite its extraordinary success, quantum mechanics has since its inception 
some seventy years ago, been plagued by conceptual difficulties. According to orthodox 
quantum theory, the complete description of a system of particles is provided by its wave 
function \l> which obeys the time-dependent Schrodinger equation. 

««|f> (1) 

According to Bohm [18], the complete description of a quantum system is 
provided by its wave function ty(q,t), q G R 3 , and its configuration Q G R 3 where Q is 
the position of the particle. The wave function, which evolves according to Schrodinger's 
equation(Eq.(l)) choreographs the motion of the particle which evolves according to the 
equation 

dQ = h 7m(tt*W) 

dt ~ m ^ ' 

where V = In the above equation H is the usual nonrelativistic Hamiltonian for 
spinless particle given as 

h 2 

H = -—V 2 + V (3) 
2m 
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Equations (1) and (2) give a complete specification of the quantum theory 
describing the behaviour of any observables or their effects of measurement. Note that 
Bohm's formulation incorporates Schrodinger's equation into a rational theory, describing 
the motion of particles, merely by adding a single equation, the guiding equation(Eq.(2)). 
In so doing it provides a precise role for the wave function in sharp contrast with its rather 
obscure status in orthodox quantum theory. The additional equation(Eq.(2)) emerges in 
an almost inevitable manner. Bell's preference is to observe that the probability current 
and the probability density p = \&*\& would classically be related by j = pv obviously 
suggests that 

^ = V = 3 /P (4) 
Bohm, in his seminal hidden-variable paper wrote the wave function \1> in the 
polar form \1> = Re %s l n where S is real and R > 0, and then rewrote the schrodingers's 
equation in terms of these new variables, obtaining a pair of coupled evolution equations, 
the continuity equation for p = R 2 as 

f = -v.(H (5) 

which suggests that p be interpreted as a probability density, and a modified Hamilton- 
Jacobi equation for S, 

— + H(VS,q) + V q = (6) 
where H = H(p,q) is the classical Hamiltonian function corresponding to Eq.(3), and 

Eq.(6) differs from the classical Hamilton- Jacobi equation only by the appearence 
of an extra term, the quantum potential V q . Similar to the classical Hamilton- Jacobi 
equation, Bohm defined the quantum particle trajectories, by indentifying VS* with mv, 

by 

dt m U 
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which is equivalent to Eq.(4). This is precisely what would have been obtained classically 
if the particles were acted upon by the force generated by quantum potential in addition 
to the usual forces. Although an interpretation in classical terms is beautufully laid down 
in the above equations, one should keep in mind that in so doing, the linear Schrodinger 
equation is transformed into a highly nonlinear equations(eqs.(5) and (6)). By taking the 
gradient on both sides of Eq.(6) we obtain 

d 1 

— v = -(v.V)v - v x (V x v) V(V + V q ) (9) 

at m 

Defining the quantum current as j(q,t) = ~Im[^*(q,t)W(q,t) - ty(q,t)W*(q,t)] = 
p(q,t)v(q,t) and using the equation V x v = we readily obtain the expression for the 
motion of the quantum current as 



-j = -w(v.j) - (j.v)v - ^v(y + v q ) (io) 

Eqs.(5), (6), (9) and (10) describe the motion of a quantum particle in the 
hydrodynamical representation of TDSE. However, the many-particle description of the 
BQH can be found elsewhere [24]. It may be noted that density alone cannot sufficiently 
describe a quantum system, one requires both density and the quantum current for the 
purpose. As is evident, the motion of a quantum particle is governed by the quantum 
current vector j unlike the TDSE where the time propagator e lHt has the key role for the 
particle's motion. The difficulties arising out of the evaluation of the exponential of an 
operator in more than one dimension is completely bypassed in the HDE. Although the 
hydrodynamical equations resemble the classical fluid dynamical equations, the quantum 
identity is prevailed because of the fact that the quantum current evolves with respect to 
a potential V q which has no classical analogue [19]. It should be be noted that the term 
V q was inherrently present in the expression to stabilize the hydrodynamical approach 
to the TDSE. The numerical instability in the hydrodynamical approach to the TDSE 
without the presence of V q term may be related to the "shock" formation in the classical 
hydrodynamics(cf. Navier Stokes equation) without some fictitious smoothing potential. 
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In the numerical solution we shall work with the equations governing the motion of den- 
sity(Eq.(5)) and the quantum current (Eq.( 10)). The motivations of considering the above 
equations lie in the fact that (i) density p and quantum current J are uniquely defined for 
a given potential in many-body system whereas the phase S can be multivalued(S=S± 
mr, n=even) and (ii) they are quantum mechanical observables. The equations 5 and 6 
suggest that one can obtain density and the quantum current directly for t > provided, 
the values were known at t — 0. Thus, the scheme by-passes the evaluation of the wave 
function during the the occurrence of the dynamics. However, at t=0, one has to solve 
the time-independent Schrodinger equation for the wave function and calculate p(q, 0) 
and j(q,0). 

3 Quantum Optimal Control and Bhomian Quantum 
Hydrodynamics 

Quantum optimal control theory seeks the design of an external field to fulfill 
a particular objective. This section will provide the rigorous mathematical formulation 
of the hydrodynamic method to design an optimal time dependent field that drives a 
quantum wave packet to a desired objective at the target time t=T. For this purpose, 
consider a general target expectation value defined as &t = Jo Qp(x, T)dx, where 9 is an 
observable operator and p(x, T) is the probability density which obeys the hydrodynamical 
equations, viz., Eqs.(5) and (10). The goal is to steer Ot as close as possible to a desired 
value Q d . We define a quadratic cost functional as 

J q = l -u a {Q T - Q d f (11) 

Minimization of J q amounts to the equalization of Ox to Q d . However, p in the above 
equation must obey the hydrodynamical equations, viz., Eqs.(5) and (10). Thus, we have 
to fulfill this constraint whereby we obtain the unconstrained cost functional as 

J = J,-ljx l{x , t)l ^ + a J^l ]dxdt (!2) 



7 



- / / A2(x ' + ^ ( 7 } + p TTx {v + v « + v ^^ dxdt 

where V ext (t) represents the external potential due to the interaction between the particle 
and the electric field, E(t) to be designed. 

Thus, in the above equations, we have introduced two Lagrange's multipliers 
Xi(x,t) and \ 2 (x,t). There exists another constraint involving the total energy in the 
field which must be imposed on the optimization procedure. This constraint takes the 
form 

\u e [( T E 2 {t)dt-E p \ =0 (13) 
l Jo 

where E p is the energy of the pulse and E(t) the field to be designed. The parameters u a 
and uj e are the positive weights balancing the significance of the two terms viz., J q and 
J e = \uj e Jq E 2 (t)dt respectively. The term J e = \u e Jq E 2 (t)dt represents the penalty 
due to the fluence of the external field. So the full unconstrained cost functional takes 
the form : 

I = Jl-//M.,t)[^ + ^ (14) 



\ 2 (x,t)[^^ + J^) + p^(V + V„ + V ext (t))]dxdt 



-II t 

+ ^e[J Q E 2 {t)dt-E p ] 



In this equation J is seen to be a functional of five functions, viz., p(x,t), j(x,t), Xi(x,t), 
\2(x,t) and E(t), all of which are real, unlike in the conventional method [6-8]. In the 
above equations capital J is for cost functional and small j is for quantum current. The 
total variation of J can be written as 

5 J = [I ^/ 5p(x, t)dxdt + / / ^/ Sj(x : t)dxdt (15) 
J J dp(x,t) J J oj(x,t) 

+ / / «3b * Al {x ' t)dxdt + 1 1 ^h) 6Hx ' t)ixdt + 1 1 m) mt)dxdt 

For any optimal solution 5 J — 0, which gives 

5J 5J 5J 5J 5J 



Sp(x,t) Sj(x,t) 5Ai(x,i) SX 2 (x,t) 5E(t) 



(16) 
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We have provided in Appendix-A the full expression for 5 J. Comparing Eq.(16) 
Eq. (A. 14) (see Appendix-A) we obtain 



5\i(x,t) dt dx 

5J =-^-^)-^+v t+ v -W ) = o 



5X 2 (x, t) dt dx p dx 



5 J d\ 2 d\i d 7. X 2 dj X 2 j dp 



5j(x,t) dt dx dx p p dx p 2 dx 

SJ dX 1 X 2 j di d 7 2 . X 2 j 2 dp 

+ 2^--^- -— {X 2 —)-2- - ' 



5p(x, t) dt p 2 dx dx p 2 p 3 dx 

+ A^h pll2 Yx {X2p) = {} 



' )J u a [Q T -Q d ]x-X 1 (x,T)=0 



5p(x,T) 

5.J 



6j(x,T) 

5 J r d 



-X 2 (x,T) = 



/o 
X 2 (x,t)p(x,t)—fj,(x)dx + u! e E(t) = 
dx 



5E(t) 

Eq.(19) and (20) can be rewritten in a simple form as 

d\ d 

-^ + -^(X 2 v x ) + S 1 [p,j,X 2 ]=0 

and 

d X d d 

~bT + d^ {XlVx) ~ x ^ v + Vq{X2) + Vext) + S2[pJ ' X2] = 

where 

a = -2^ 

p dx 
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and 



S 2 = - 



X 2 ^-(V q (p)-V q (X 2 )) 



J 



fd\2 

? 2 dx 



(27) 




Note that the above expression for O-r restricts the operator G being only 



a multiplicative operator, for example, the distant x which we have used in the sub- 
sequent numerical calculations. However, other forms of operator can also be consid- 
ered in the BQH-QOC formulation with the different constraint expressions, e.g., if is 
the momentum operator(j>) we would require the constraint equations (5) and (9) since 



being the extra source terms Si and S 2 . The source terms depend on p and j. v \ in 
the above equations is the velocity associated with the Lagrange's multiplier and is given 
as v x = ^ and V q (X 2 ) is given by V q (X 2 ) = Notice that Eqs.(17) and (18) 

are the equations of motion for the probability density and the quantum current density 
respectively obtained in section 2. Whereas Eqs.(24) and (25) are the equations of motion 
for the two Lagrange's multipliers A2 and Ai respectively. It should be noted that in 
obtaining the above equations(see Appendix-A) we have assumed no variation on either 
p(x, 0) or j(x,0). Thus, we start from an initial(t=0) p(x, 0) and j(x,0) to solve Eqs.(17) 
and (18) for p(x,t) and j(x,t) respectively. Equations (24) and (25) can be solved for 
\i(x,t) and X 2 (x, t) provided a starting value Xi(x, t s ) and X 2 (x, t s ) were known. These 
have been obtained from Eqs.(21) and (22) respectively as 



where t s = T, the final time. Thus, one has to perform backward propagation for solving 
both the equations of motion involving Xi(x, t) and X 2 (x, t). Having calculated p(x,t), 
j(x,t), Xi(x, t) and X 2 {x, t) as described above, one has to carry out an optimization of the 



Pt = m Jo p(x,T)V S(x,T)dx. 

The equations for Ai and A2 ressemble to that of p and j with only difference 



Xi(x, t s ) = 0J a [Qx — ® d ]x and X 2 (x,t s ) = 



(28) 
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quadratic cost functional(Eq.(ll)) with respect to the electric field E(t) which, according 
to Eq.(23), takes the form 



This constitute the details of the BQH-QOC method. 

4 Application to HBr Molecule 

We have said in the preceeding section that we needed the initial density p(x, 0) 
and the quantum current j(x,0) in the presnt method. These have been evaluated by 
solving the time independent Schroedinger equation for HBr molecule in the 1 J2 + state 
where the the potential energy is assumed Morse type of the form [20] 



where (3 = cj e ( 1 ^) 1 / 2 , D e = with u e = 2648.975cm" 1 , u e x e = 45.217cm" 1 , x e = 

1.41443 angstrom and \x being the reduced mass of HBr. 

Having obtained p(x, 0) and j(x,0) we carry out the control by the present method. 
Followings are the necessary steps for the computer implementation of the present method 

A. Present Method : 

Step 1: Make an initial guess for the electric field E(t), which is zero in our 
calculation. 

Step 2: Solve the coupled equations, viz., Eq.(17) and (18) for p(x,t) and 




(29) 



V = D e (l - exp{-(3{x - x e ))) 2 



(30) 



j(x,t) respectively starting from p(x,0) and j(x,0). The solution is 
done by using the Flux-corrected transport (FCT) algorithm [21] mod- 
ified by us for the purpose of solving the quantum hydrodynamical 
equations [16,17]. In doing so, we adopt the Eulerian scheme 



Step 3: Evaluate the final values for X±(x, T) and X 2 (x,T) given by Eq.(28). 
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Step 4: Use Ai(x,T) and A 2 (x,T) for solving Eqs.(24) and (25) for Ai(x,t) 
and A2OE, t) respectively. This is done by backward propagation, by 
putting dt=-dt(see ref.16). We follow the same method as in step 
2 for solving these equations. It should be noted that Eqs.(24) and 
(25) have source terms which depend on p(x,t) and j(x,t) calculated 
in step 2. 

Step 5: Calculate the quadratic cost functional given by Eq.(ll). 

Step 6: Optimize the function in Eq.(ll) with respect to the electric held, 

E(t) given by Eq. (29). Here we use the conjugate direction search 

method [22] for the optimization. 

Step 7: Iterate step 2 to step 6 until a convergence criterion is satished. 

The external potential is of the form V ext (x,t) = —/i(x)E(t), where /i(x) is 
the dipole function for HBr and is given by [23] /i(x) = A + A\(x — x e ) + A 2 (x — x e ) 2 
where A = 0.788, A x = 0.315 and A 2 = 0.575. In our calculation the range of spatial 
dimension is < x < 12 a.u., that of time is < t < 2000 a.u. Total number of spatial 
mesh points is 60 which gives Ax = 0.2 a.u. Similarly, total number of time steps is 2000, 
which corresponds to At = 1.0 a.u. u e in Eq.(27) is taken as 0.5, and oo a as 1000. The 
target operator is 9 = x and Q d = 3.0 a.u. and 3.5 a.u. 

Figure 1 shows the electric fields corresponding to two different values of 
Q d viz., 3.0(solid lines) and 3.5(dotted lines). These pulses excite several vibrational 
states(not shown here) mainly by a sequence of single quantum transitions. The peak 
value of the field is ~ 0. 08a. u. (corresponding intensity is ~ 10 14 Wcm~ 2 ) for Q d = 3.5a.u. 
and ~ 0. 02a. u. (corresponding intensity is ~ 10 13 Wcm~ 2 ) for Q d = 3.0a.w.. The detail 
characterization of the optimal field can however, be made by Fourier transforming the 
field. Fig.2 shows the average distance < x > as a function of time. Notice the desired 
control of < x >= 3.0 and 3.5 a.u. at T=2000 a.u. is obtained through the oscillatory 
motion of the packed induced by the optimal electric pulse(Fig.l). Figure 3 shows the 
initial and final densities for the two values of Q d . The packet is distorted in shape while 
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approximately retaining its original variance during the evolution. During the optimiza- 
tion process the total integrated probability density remained at unity up to a deviation 
of 10~ 7 . The number of iterations in the optimization to achieve the results is 5 and it 
takes only 3 minutes(real) on a IRIX IP30 machine with R4400 6.0 CPU. As a test for the 
acceptability of the present method we have carried on the following experiment : The 
electric fields(Fig.l) so obtained have been pluged into the TDSE and then solved for the 
wave function. The results for the density and the expectation value of < x > resemble 
accurately to that given in Fig. 2. 

5. Conclusion 

In the present paper we have presented a new scheme for carrying out the optimal 
design based on BQH. We have derived the control equations to obtain a time dependent 
external field with an illustration for the manipulation of the vibrational motion of HBr 
molecule in the 1 J2 + state. The working dynamical variables in the BQH , viz., p(Fig.3), 
j, Ai and A2 are relatively slowly varying spatial functions(Fig. 4) compared to the wave 
function(fig.4, curve a) which apparently enhances the efficiency and the numerical saving 
of the BQH-QOC method for controlling dynamics. 

Although the illustration of our new method has been made in one spatial 
dimension, the approach is general and is directly extendable to higher dimensions and 
a wave packet dynamics in four dynamics has already been performed [16] within our 
method. The use of the alternating direction implicit (ADI) [16,17] in the present method 
makes the quantum control calculation much easier compared to the conventional method, 
especially for the multidimensional problem. In the conventional optimal control theory, 
the role of the complex Lagrange's multiplier is to provide feedback[6] for designing the 
electric field and guide the dynamics to an acceptable solution. The BQH-QOC method, 
on the other hand, introduces two such Lagrange's multipliers, Ai and A2 both of which 
are real variables. The first Lagrange's multiplier Ai, which corresponds to the quantum 
current j(cf. Eq.(18) and (25)) has however, no direct role to provide feedback for designing 
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the electric field(Eq.(29)) and only guides the dynamics in conjunction with the second 
Lagrange's multiplier A 2 . It may be worth mentioning that since the quadratic cost 
functional(Eq.(ll)) is a functional of density, the Lagrange's multiplier \ 2 (equivalent to 
the density p, cf. Eq.(17) and (24)) enters into the expression for the optimal electric 
field(Eq. (29)). However, cases where one desires to manipulate the quantum flux(which 
is directly related to the quantum current j) by constructing a quadratic cost functional 
dependent on j, the Lagrange's multiplier Ai will appear explicitly into the expression for 
the external field. 



where the dynamics may lead to the creation of the nodes in the density profile since the 
quantum potential appearing in the constraint equation blows up in the occurence of such 
an event. However, such occurence of nodes can be countered by fixing a lower limit to the 
density of the order of the machine precision. This in other words means that one never 
encounters an absolute nodal point where the density is exactly zero. Future studies 
need to explore the other area of control within the BQH-QOC method, for example, 
controlling the quantum flux. 
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The variation of J given by Eq.(14) has to be taken with respect to p(x, t), j(x,t), Ai(x, t), 
\2(x,t) and E(t). Any variation 5p(x,t), 5j(x,t), 5Xi(x,t), 5\2(x,t) and 5E{t) will lead 
to the variation 5 J given as 



It should be pointed out that the present method could prove hard in cases 



APPENDIX 




dp dj 



(Al) 
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- / / x 4 Sj + s<7* " 7' 5p) + sp rJ v + 1 « + v -> 

+p— (SV + SV q + 5V ext )\dxdt + u e [ E(t)6E(t)dt 
ox J 

Now, we have SV = and SV ext can be written as 8V ext = —5(p(x)E(t)) = —fi(x)5E(t) — 
E(t)S/i(x). Since /i(x) is kept fixed, we get SV ext = —/i(x)SE(t). J q in the above equation 
is given by 

J q = \w a [ J p(x,T)xdx - x d cm ] 2 (A.2) 

Hence 

5J g = uj a \<x> (T) - x d cm ] J x5p(x, T)dx (A3) 
Substituting Eqs.(A.3) into Eq.(A.l) we obtain 
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SJ = w„[< z > (T) - xf j j xSp(x, T)dx - j J[^+ ^SMdxdt 

-//'I + S ( 7 )+p S (v ' + ^ + ^» A 

— J J \i—5pdxdt — j j X\—5jdxdt 

- j j X 2 -^5jdxdt-2 J J X 2 --^5jdxdt (A4) 

-2 / / X 2 -^-5jdxdt + 2 f [ X 2 ^%b 2 dxdt 
J J pox J J p 2 ox 

J J ^ 2 ^2~^3pdxdt + 2 J j X 2 ~^bpdxdt 



+ 
-2 



/ l X2 ^~L 5f)dXdt ~ l J X ^( V + V 9 + V ext)Spdxdt 

- j J X 2 p-^5V q dxdt + J J X 2 p-^(p(x)5E(t))dxdt 
+u e J E(t)5E(t)dt 

The 4-th and 6-th terms in the above equation can be simplified by integration by parts 
as follows 

I J Xl ~Q£^Pd x dt = J Ai (x, T)8p(x, T)dx — J X±(x,0)5p(x,0)dx 

- j j ^5p(x,t)dxdt (A5) 



d 

djdxdt = I X 2 {x,'l ) 

r r r>\,, 

-Sj(x,t)dxdt (A6) 



J j X 2 —5jdxdt = J X 2 (x,T)5j(x,T)dx — J X 2 (x,0)5j(x,0)dx 

dX 2 „ 



dt 

Terms 5-th, 7-th and 10-th can similarly be expressed by the integration by bparts as 
follows 

/ J ~oh^ dxdt = J ^( x ^ t )^j( x r,t)dt - j X 1 (x h t)Sj(x h t)dt 
-J J ^-5j(x,t)dxdt (A.7) 
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A 2 --^-5jdxdt 
p ox 



X 2 (x r ,t)j(x r ,t) 



p(x r ,t) 

SSw^ dxdt 



Sj(x r , t)dt 



\ 2 (xi,t)j(xi,t) 
(A8) 



5j(xi,t)dt 



J 2 d 
X 2 ——5pdxdt 
p A ox 



\ 2 (x r ,t)j 2 (x r ,t) 

p 2 (x r ,t) 

J J T x (x 4 )S » dxdt 



5p(x r , t)dt — J 



^2{xut)j 2 {x h t) 
(A9) 



Sp(xi, t)dt 



Term 15-th is 



/ / X ^P^. ( y^ x ) 5E( y t )) dxdt = j J Mp-7^p(x)5E(t)dxdt 



d 



(A10) 



14-th term involves the variation in J due to the change in the quantum potential 5V q , 

n 2 vV /2 



where V q is given by V q = — ^ f /2 . This gives 



h 2 d 2 



1 



5Vq App^dx 2 v p 1/2 



Sp) 



1 



d 2 



(All) 



4ppV 2 dx 21 

By the integration by parts we simplify the 14-th term as follows 
/ J X2 P~^dV q dxdt = J \ 2 {x r ,t)p(x r ,t)5V q (x r ,t)dt - J \ 2 (xi,t)p(xi,t)5V q (xi,t)dt 



+ 



+ 



d . d . . 1 . 1 . . 

( — (A 2 p) 77777^ ) 77777^ |x, <W« , t) dt 



(A12) 



dx K dx K "^'2p 1 l 2 '2p 1 / 2 ^ 1 

Where x r and :r; are the right and left ends of the one dimensional grid, and 
F(x r ) — F(xi) where F(x) is any function. The first and the second terms in 



Eq.(A.12) are the contributions due to the change in the quantum potential at the two 
ends of the boundary only. Since, we take a large grid, p at the two ends of the grid 
are very small and can be assumed constant. This leads to V q (x r ,t) and V q (x h i) being 
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very high constant values at any time and hence 5V q (x r ,t) = SV q (xi,t) = 0. With the 
same argument we can also neglect the contributions due to the terms 3rd, 4th and 5th. 
Combining all the terms we obtain the full variation in J as 



5 J = co a [< x > (T) — x d cm ] J x5p(x, T)dx -J J[^ + ^SXidxdt 

-J X 1 (x,T)5p(x,T)dx + J X 1 (x,0)5p(x,0)dx (A13) 
J J ~£^~^Pd x dt ~ j X 2 (x, T)5j(x, T)dx + J X 2 (x, 0)5j(x, 0)dx 
J f ~^Sjdxdt- j \-i(x r ,t)8j(x r ,t)dt + J X 1 (xi,t)5j(xi,t)dt 



+ 
+ 



+ // ^5j(x,t)dxdt-2 J ^^^ 5j(x r ,t)dt 

+2 / X ^ Xu f )j ^ ht) 5j{x h t)dt+ [ [ -^-(X 2 3 -)5j(x,t)dxddt 
J p{xi,t) J J ox p 

f X 2 (x r ,t)j 2 (x r ,t) f X 2 (x r ,t)j 2 (x r ,t) 

+ J 7{x-T) 6 P(^t)dt-J 7( - J) 5p(x r ,t)dt 

-J J ^(X 2 J -)8pdxdt-2 J J \ 2 -^-8jdxdt 

+2 f f \ 2 L^-8jdxdt+ f f X 2 %%bpdxdt 
J J p l ox J J p 2 ox 

-2 J J \ 2 3 1^5pdxdt- J J \ 2 ^-(V + V q + V ext )Spdxdt 

1 r r d 2 . d . 1, 1 - , , 

J J w^^w^W"* ' 
+ -J ! ii M w^ pinSfixdt 

+ J J X 2 p-^p(x)5E(t)dxdt 

This expression has 26 terms. Out of which, 5-th and 8-th terms can be dropped 
because we do not vary the initial density and quantum current. Again, 10-th, 11-th, 
13-th, 14-th, 16-th and 17-th terms can also be dropped with the assumption that p(x,t) 
and j(x,t) are very small at the boundary. Thus, the actual full variation in J becomes 
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5J = uj a [< x > (T) — x d cm ] J x5p(x, T)dx - j J[^ + ^SXtdxdt 



dxdt 



- J X 1 (x,T)5p(x,T)dx + j j^-5pdxdt (A 14) 

- J X 2 (x,T)5j(x,T)dx + J J ^5jdxdt 

+ J j *^-5j(x,t)dxdt + J J -^-(X 2 —)Sj(x,t)dxddt 

-f J -^{X 2 ^)5pdxdt-2 J J X 2 ~-5jdxdt 
+2JJ X 2 j^^5jdxdt + j j X 2 ^^5pdxdt 



-2 J! X 2 ^^6pdxdt- 1 1 X 2 -^-(V + V q + V ext )Spdxdt 



1 r r d 2 d 1 1 

'J J W [ ^ M ^ ] 2p^ Spdxdt 

1 r r d x 1 d 2 , , / / , d 



+ pll Tx {X2p) Ap^ 2 W pl,25pdxdt + II ^P^(x)6E(t)dxdt 
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Figure Captions 

Figure 1: Optimal electric Geld shown as a function of time for Q T = S.Oa.u. (solid 
line) and Qt = 3. 5a. u. (dotted line). 

Figure 2: The expectation values < x > shown as a function of time for Qt = 
3.0a.w. (solid line) and Qt = 3. 5a. u. (dotted line). 

Figure 3: lnitial(t=0) (dotted line) and hnal(t=T) (solid line) density correspond- 
ing to Qt = S.0a.u.(lebel a)) and Qt = 3.5a.u. (label b)). 

Figure 4: Hydrodynamical variables, viz., j(x,T)(b), Xi(x,T), \ 2 (x,T)(c) and 
the real(solid) and imaginary(dotted) values of the wave function(a) 
plotted as a function of x. Notice that the hydrodynamical variables 
are smooth spatial function unlike the wave function. 
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